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Abstract—This paper is devoted to the application of the 
Semi-Analytical Cloud Retrieval Algorithm (SACURA) to the 
cloud-top height determination using data from the Global Ozone 
Measurement Experiment (GOME) instrument onboard Earth 
Remote Sensing satellite (ERS-2). In particular, measurements 
of the top-of-atmosphere reflectance in the oxygen absorption 
A-band are used. The technique is based on the asymptotic 
radiative transfer theory as applied to studies of the oxygen 
absorption bands in reflected light. Our approach is valid for 
optically thick clouds with cloud optical thickness larger than 
approximately 5. The accuracy of the algorithm is checked against 
independent retrieval techniques for completely cloudy pixels. 
In particular, the close coincidence with data obtained from the 
Along Track Scanning Radiometer (ATSR-2) onboard ERS-2 is 
found. Some results of retrievals using these different instruments 
disagree (up to 2 km). This is explained by the strong horizontal 
inhomogeneity of clouds under investigation, which is not ac- 
counted by the SACURA or, possibly, by well-known problems of 
infrared techniques as applied to low-level clouds. The effective 
cloud geometrical thickness | is also retrieved. This parameter is 
defined as the geometrical thickness of a vertically homogeneous 
cloud, which allows for the minimization of differences between 
observed and calculated top-of-atmosphere reflectance spectra. 
For inhomogeneous clouds, the value of | differs from a real cloud 
geometrical thickness, but it gives us an indication of the possible 
existence of the multilayered cloud system in the field of view of 
the optical instrument. 


Index Terms—Clouds, radiative transfer, remote sensing, spec- 
troscopy. 


I. INTRODUCTION 


OST OF THE modern satellite-based spectrometers 

and radiometers have a capability to get information on 
a cloud-top height distribution on a global scale [6]-[10], [14], 
[17]-[19], [26]. As far as measurement techniques and retrieval 
methods are concerned, they differ considerably both in terms 
of their theoretical basis and accuracy, depending on the instru- 
ment. They range from active measurements, using lidars and 
radars on satellite platforms [20], to passive techniques, based 
on the processing of information, contained in the thermal 
infrared radiances [7] or reflected light polarization [8]. Stereo 
photogrammetric methods [16] and methods based on the 
measurements of the top-of-atmosphere (TOA) reflectance in 
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the oxygen A-band [10], [17] are also very popular techniques 
nowadays. 

The subject of this paper is the first application of the 
improved oxygen A-band-based cloud-top height retrieval 
technique developed by Rozanov and Kokhanovsky [24] to the 
Global Ozone Monitoring Experiment (GOME) measurements. 

The physics behind the technique is quite simple. Indeed, let 
us assume that we have a gas in a planetary atmosphere, which 
absorbs almost all incident radiation in a narrow spectral band. 
Then, the difference in the reflection between wing and center 
of this band (depth of a band), measured by a receiver on a 
satellite orbiting the planet, will depend on the cloud altitude. 
Gas concentrations generally decrease with the distance from 
the ground. Clouds at a high altitude do not allow photons to 
penetrate to low atmospheric layers and be absorbed there. So, 
the depth of a spectral band will decrease, if clouds are present 
in the field of view of a sensor. 

However, the absorption of photons takes place not only in 
the atmosphere above the cloud but in the cloud itself, too. This 
leads to the enhancement of the absorption. Thus, the depth of 
the absorption band depends both on the cloud altitude and the 
cloud geometrical thickness. 

Indeed, the very first results of the cloud-top height retrieval 
from the measured reflectance spectrum in Oz A-band [25] 
have shown that this enhancement of the absorption has to be 
taken into account to derive an accurate value of the cloud-top 
height. 

Therefore, the proposed Semi-Analytical Cloud Retrieval 
Algorithm (SACURA) allows for the retrieval of both the 
cloud-top height h and the effective cloud geometrical thick- 
ness l simultaneously. This is the main new feature of the 
SACURA as compared to other cloud retrieval techniques. 
Note that the retrieved value of | corresponds to the case of a 
homogeneous cloud layer. Therefore, it can be considerably 
biased as compared to a sum of the geometrical thicknesses 
of cloud layers in multilayered cloud systems. However, the 
account for the effective value of J increases the accuracy of 
the cloud-top height retrieval (see below). 

The solution of the inverse problem is very complicated in the 
case studied. Popular lookup-table techniques are not applicable 
due to a large number of parameters involved. The minimization 
techniques based on the exact solution are very time consuming. 
They do allow for single runs but cannot be considered as an 
option for operational retrievals, due to the computation speed 
constrains. 

For that reason, our retrieval algorithm uses the approximate 
analytical solution of the radiative transfer equation given by 
Kokhanovsky and Rozanov [11] to speed up the retrieval con- 
siderably. This is another new feature of the SACURA. 
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As it follows from our theoretical studies [24] the accuracy 
of the SACURA is quite high for homogeneous cloud layers. It 
was also established that the vertical inhomogeneity of clouds 
such as vertical profile of liquid water content and especially 
multilayered cloud systems can influence the accuracy of the re- 
trieval considerably. The maximal errors of the determination of 
the cloud-top height occur if the observed cloud field is remote 
from the homogeneous assumption made in the forward model 
(e.g., the presence of multilayered clouds). The retrieved value 
of J can be in such cases so large that a homogeneous cloud 
must have a value of the geometrical thickness larger than h. 
This points out to the fact that the retrieved geometrical thick- 
ness can be used as a quantitative criterion to judge whether our 
homogeneous assumption corresponds to the real atmospheric 
situation. 

Note that the retrieved cloud-top height does not give any 
indication on the correspondence of the forward model to a real 
atmospheric situation. Most modern retrieval techniques have 
this weak point. 

The performance of our semianalytical cloud retrieval algo- 
rithm is highly dependent on the characteristics of the instru- 
ment used. In particular, high spectral and spatial resolution 
in oxygen A-band is needed. Unfortunately, such data are not 
available to us. An instrument [the Orbiting Carbon Observatory 
(OCO)] to fully explore possibilities of our technique has been 
already designed [2], [13]. It has 1.0 x 1.5 km? spatial resolu- 
tion. The resolving power of the OCO is equal approximately to 
17500. It translates into approximately 0.04 nm for the spectral 
resolution. However, the OCO was not launched in space so far. 
So, we used TOA reflectance measurements from the GOME in- 
strument [1] onboard Earth Remote Sensing satellite (ERS-2). 
The spectral resolution of this instrument suits our needs. The 
GOME instrument has the resolution ~0.33 mm in the oxygen 
A-band. This allows to obtain approximately 60 measurement 
points in the spectral range 758—771 mm. The spatial resolution 
is, however, very poor (typically, 40 x 320 km?). It means that 
the selection of a pixel completely covered by clouds is not a 
very easy task. Indeed, we need to avoid broken cloud cases. 
Fortunately, the GOME makes measurements at the enhanced 
resolution mode (40x 80 km”) 3 times per month. So, we mostly 
used these rare occasions to see the performance of our tech- 
nique as applied to experimental data. 

We have selected a single GOME pixel with complete cloudi- 
ness at 43.23°N, 27.44°W at 12:52 UTC on October 15, 1996 
above the Atlantic ocean. The cloud-top height for this pixel was 
estimated from other measurements as well (see details in [4]). 
The retrieval for this case is considered in the Section II. 

Section II is devoted to the comparison of GOME retrievals 
with independent data obtained from the Along Track Scanning 
Radiometer-2 (ATSR-2) based on infrared measurements. The 
ATSR-2 and GOME instruments are on the same space plat- 
form. However, they have different spatial resolution and dif- 
ferent orientation of footprints. In particular, ATSR scans along 
a satellite track and has a capability of dual viewing with the 
spatial resolution 1 km?. So, the comparison should be consid- 
ered as a qualitative one. However, it does provide us with the 
valuable information. This is due to the fact that the cloud-top 
altitude usually does not change very much for extended cloud 
fields. 
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Fig. 1. Normalized spectral reflection function R(X) in O2 A-band measured 


by the GOME instrument (symbols). The size of the pixel is 40 x 320 km?. 
Measurements have been performed at 12:52 UTC on October 15, 1996 above 
the Atlantic ocean (43.23°N, 27.44°W). The solar and observation angles are 
equal to 53° and 0°, respectively. Solid line denotes calculated results. 


II. APPLICATION OF THE SACURA TO GOME DATA: 
A SINGLE-PIXEL STUDY 


In this section we use a simplified version of the SACURA 
based on the one-parameter retrieval algorithm (see details in 
Appendix). This allows us to show the dependence of the re- 
trieved cloud-top height on the geometrical thickness and to 
demonstrate that the geometrical thickness can be evaluated as 
an additional parameter. 


A. Determination of the Cloud-Top Height Using Experimental 
Reflectance Spectra 


Let us derive the cloud-top altitude from the reflectance 
spectrum R(A) measured by the GOME in the spectral interval 
758-771 nm (O2A-band) over Atlantic Ocean. The measured 
reflectance spectrum normalized by the value R* = R(758 nm) 
is shown in Fig. 1 by symbols. We see that the absorption in 
A-band is much more stronger than outside of the band . As a 
matter of fact, the reflectance in the center of the oxygen band 
for a clear sky case is almost zero. It is around 0.2R* in the 
case studied. This gives us a clear indication that a cloud field 
did present in the field of view of the instrument during the 
time of measurements. 

The same GOME cloudy pixel has been already studied in 
respect to the cloud-top height determination by de Beek et al. 
[4]. However, they used another spectral region (392-395 nm) 
and a different retrieval technique (Ring effect). They report that 
the cloud top has been located at 7 + 1 km for the case studied. 
De Beek et al. [4] confirmed this finding using independent in- 
frared measurements by Meteosat, recorded at 12:00 UTC on 
October 15, 1996 with account for the temperature measure- 
ments of a balloon radiosonde, launched by the research ship 
Polarstern at 42°N and 35°W at 10:18 UTC on the same day. 
So, we are pretty sure in the value of h for the case studied. 

According to the GOME cloud coverage data, the pixel 
studied was completely covered by clouds. Therefore, the 
SACURA should be applicable in this case. Before starting 
the retrieval algorithm, we need to define some cloud and 
atmospheric parameters. 

So, we use the value of the cloud optical thickness 7 obtained 
from the measurements in the visible (r = 27.7) in the assump- 
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geometrical thickness from the measured GOME spectra given in Fig. 1. The 
extrapolation of the results to / = 0 is shown with dotted line. 


tion that the vertical profile of the liquid water content (LWC) 
is constant. Also here and below we assume that ice particles do 
not present in clouds studied and the droplet size distribution is 
given by the Cloud C1 distribution [5] with the effective radius 
of droplets equal to 6 um. The atmospheric aerosol model co- 
incides with that used by Kokhanovsky and Rozanov [11]. 

The cloud geometrical thickness / is not known a priori. 
Therefore, we estimated the value of h for a number of discrete 
values of / ranging from 1-8 km with the step 1 km. The depen- 
dence of the retrieved value of h on the geometrical thickness 
l is shown in Fig. 2. We see that this dependence is well 
pronounced and close to the linear function. It follows from 
Fig. 2 that the retrieved value of h increases with J. This can 
be easily understood. Indeed, larger values of the geometrical 
thickness of a cloud lead to larger light absorption by a cloud. 
It means that the increase in the value of J leads to the more 
pronounced minimum of the modeled reflectance spectrum. 
The depth of the band is fixed, however. It is obtained from 
measurements. Therefore, to fit measurements for larger values 
of l, the algorithm produces larger values of h and, therefore, 
less absorption above cloud. This effect is clearly seen in 
Fig. 2. The obtained result shows that for accurate retrieval of 
the cloud-top height we need to estimate the effective cloud 
geometrical thickness. 

It is reasonable to take as such an estimation the value of J, 
which provides the best correspondence of measured and mod- 
eled spectra. The cost function (see the Appendix) is quanti- 
tative measure of such a correspondence. Therefore, we must 
choose a pair (h,/), which provides the minimum of the cost 
function. 

The dependence of the cost function on the geometrical thick- 
ness is shown in the Fig. 3. We see that the cost function reaches 
a minimum for / equals to 7 km. Using data shown in Fig. 2, we 
find that the corresponding value of h is equal to 7.2 km. So, the 
solution is the pair (7.2 km, 7 km) for the case studied. 

Fig. 1 gives the comparison of the measured and calculated 
spectra in the oxygen A band. Results of calculations correspond 
to h = 7.2 km and / = 7.0 km found from the retrieval tech- 
nique. Close coincidence of calculations and measurements is 
found. Therefore, the application of the radiative transfer equa- 
tion with these values of pair (h, l) (together with derived value 
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of 7 = 27.7) can explain the experimentally measured spectrum 
quite well. 

The obtained result is worth to be discussed in more details. 
Using the retrieved value of / and the value of the cloud optical 
thickness given above, we can estimate the mean cloud extinc- 
tion coefficient. The cloud extinction coefficient defined as the 
ratio T/L is equal to ~ 4 km" for this case. This value is un- 
realistically low. This leads us to the conclusion that a single 
homogeneous layer model used in the retrieval is not represen- 
tative for the case studied. 

Although the retrieved value of the cloud-top height agrees 
well with independent retrievals from other techniques, this 
could be a result of a pure coincidence. Therefore, more 
retrievals should be done to establish the technique presented 
here. Then, independent measurements of h should be also 
performed. Unfortunately, such measurements (e.g., using laser 
systems) are rare. 

Figs. 1-3 justify our algorithm for the cloud-top height deter- 
mination. However, the problem of uniqueness of the solution 
arises. Indeed there is a probability that such a good coincidence 
of satellite measurements with the data obtained from the radia- 
tive transfer theory as shown in Fig. 1 can be achieved assuming 
that the cloud layer is not uniform in the vertical direction (but, 
probably, for other values of h, l). Actually, for most of cases of 
extended cloudiness one or more cloud layers may exist. 


B. Retrieval of Cloud Top-Height Using Synthetic Reflectance 
Spectra for a Two-Layered Cloud System 


To check the assumption on the multilayer structure of the 
cloudiness in the observed scene, we have decided to perform a 
numerical experiment. To model the situation, we have used the 
synthetic reflection function obtained for a two-layered cloud 
system (in the assumption that both layers (see Fig. 4) are com- 
posed of water droplets with the constant liquid water content) 
in our retrieval algorithm for a single homogeneous cloud case. 
Then, the technique based on the A oxygen band was applied. 
Note that the upper layer is located at heights 5—8 km, and the 
lower level cloud extends from 2 until about 3 km above the 
ground surface. 

Fig. 4 shows the retrieved values of h for an effective homo- 
geneous cloud for geometrical thicknesses in the range 1-7 km. 
The retrieved value of 7 coincides with the sum of values of 7 
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geometrical thickness from simulated data for a two-layered cloud system. 
The optical thickness of the upper cloud is 7. The optical thickness of the 
lower cloud is 30. The solar and observation angles are equal to 53° and 0°, 
respectively. Upper cloud is positioned between 5 and 8 km. The lower cloud 
layer is located between 2 and 3 km. The extrapolation of the results to ł = 0 
is shown with dotted line. 
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of upper and lower layer (7 = 37). It follows from Fig. 4 that 
the dependence h(l) closely follows that for the case given in 
Fig. 2, where experimental spectra are treated. 

The dependence of the cost function on the value of l is given 
in Fig. 5. The minimum of the cost function is reached at | = 
6 km. The depth of the minimum of the cost function is much 
more pronounced than it was in the case of experimental spectra 
(see Fig. 3). The reason is that the synthetic spectra have been 
used in the case given in Fig. 5. 

Using data shown in Fig. 4 we find that corresponding value 
of h equals to ~6 km for] = 6 km. We see that the retrieved uni- 
form cloud has the geometrical thickness equal h. Taking into 
account that 7 = 37 we find that the effective cloud extinction 
coefficient is ~ 6 km. This is similar to our finding for experi- 
mental data, where also values of h and l were almost identical. 
Also we had too low values of cloud extinction coefficient in 
that case. 

Comparing obtained cloud-top height with the “true” value 
(see Fig. 4), we see that the retrieved value of h is lower by 
2 km. 

All modern cloud retrieval algorithms based on the O2A 
absorption band ignore the influence of the cloud geometrical 
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thickness on the retrieval. Our results show that this leads to the 
bias in the retrieved values of the cloud-top height. In particular, 
it follows from Fig. 4 that h(l = 0) = 3.4 km (see dotted line 
in Fig. 4). Therefore, neglecting the value of l, we obtain the 
value of h, which is by 4.6 km lower than the actual value of 
h = 8 km (see Fig. 4). So, using l as an additional retrieval 
parameter, we can estimate the cloud-top height more than two 
times better than in in the case if we neglect the absorption of 
light inside cloud. 

Similar extrapolation of the dependence h(l) given in Fig. 2 
to 1 = 0 gives the value of h equal to 4.1 (instead of 7.1, as it 
follows from infrared Meteosat measurements, Ring effect, and 
our retrievals with account for the value of /). 

The numerical experiment performed supports our idea that 
a multilayered cloud system existed in the moment of measure- 
ments given in Fig. 1. Note that the obtained result is in a good 
agreement with the data of retrieval of h and | for two-layered 
cloudy media, which is given in the theoretical study performed 
by Rozanov and Kokhanovsky [24]. However, in the last case, 
solar zenith angles and cloud layers optical thickness were dif- 
ferent from the case studied here. 

The large value of the retrieved cloud geometrical thickness 
for two-layered cloudy media can be explained as follows. The 
absorption of photons in a two-layer system occurs not only in 
the upper and lower clouds but also in the space between them. 
This is due to the multiple light scattering and radiative transfer 
from one cloud layer to another. The total light absorptance of 
the whole system is clearly larger than the sum of the light ab- 
sorptances of both layers. The value of / in our model is a kind 
of a regulator of the absorptance level in a cloud. Therefore, 
the overestimation of l (e.g., l coincides with the distance from 
the cloud top to the ground surface or even to the level which 
is below this surface) could be a result of the presence of multi- 
layered clouds in the scene under study. The retrieval technique, 
however, is applied to the case of a single homogeneous cloud 
layer. Therefore, unrealistically large values of l retrieved may 
indicate that the real situation in the atmosphere does not re- 
semble the assumptions we made in the retrieval scheme. The 
retrieval of large values of / (e.g., larger than A) is an indicator 
of a multilayered cloud system. 

The discussion presented above show that measurements in 
the oxygen A-band can reveal information on the vertical dis- 
tribution of water in cloud systems. This possibility is worth of 
exploration in future research. 


III. COMPARISON OF CLOUD-TOP ALTITUDE RETRIEVALS 
OBTAINED USING THE OXYGEN A-BAND TOA REFLECTANCE 
MEASURED BY THE GOME AND THE INFRARED 
MEASUREMENTS OF THE ATSR-2 


We were concerned with the retrieval of the cloud-top height 
for a single (but previously well-studied [4]) cloudy scene in 
the previous section. Now we move to the application of our 
algorithm to selected cloudy pixels from several GOME or- 
bits. The results of our retrievals are compared with indepen- 
dent data from ATSR-2 instrument located at the same satel- 
lite as the GOME spectrometer. The brightness temperatures of 
cloudy pixels as measured by ATSR-2 are routinely used to ob- 
tain cloud-top pressures (and heights). The European Center for 


ROZANOV et al.: DETERMINATION OF CLOUD ALTITUDES 


8 — eo 


ATSR-2 
F x—*—x SACURA & 4 
FRESCO ~ 


F 6 46 & 


6F o o o 


Cloud top height, km 
> 
T 


H N 4 
L A 
° 
F rs , 
° 
2} fi ad 


45 50 55 60 65 70 
Latitude, degrees 


Fig. 6. Cloud-top height retrieved from GOME (FRESCO, SACURA) and 
ATSR-2 data for selected pixels from orbits 50723 095 and 50723 114 of the 
ERS-2 (July 23, 1995), having cloud cover close to 1. 


15F 4-7 
a 

F a a a ATSR-2 4 

= L xx SACURA d 

x L 4 

Z 10+ 4 

Rs L 4 
v 

£ L 4 
a 

2 f; 4 

z L 4 

E 5- m 

O F 4 


-—40 —20 (0) 
Latitude, degrees 
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having cloud cover close to 1. The pixel positions are shown on the map. 


Medium-Range Weather Forecasts (ECMWF) model (ECMWF 
1995, the description of the ECMF/WCRP Level I-A Global 
data Archive) is applied to increase the accuracy of such a con- 
version. The coincidence of retrievals confirm that both methods 
are at least consistent. The disagreement shows that at least one 
of the methods produces wrong results. 

Selected results of our retrievals as compared to the ATSR-2 
data are given in Figs. 6 and 7 for scenes located both over land 
and ocean. We assumed in the retrieval that the albedo of ocean 
is equal to 0.05 and that of land is equal to 0.3. Although the as- 
sumed albedo of the land can be biased as compared to a “true” 
surface albedo, the influence of the surface on the top-of-atmos- 
phere reflectance is low for pixels studied (cloud optical thick- 
ness is larger than 35). 

Note that the spatial resolution of the ATSR-2 differs from 
that of GOME. Therefore, we averaged the ATSR-2 data to bring 
them in correspondence with the spatial resolution of GOME. 
Only enhanced resolution (40 x 80 km?) GOME pixels com- 
pletely covered by clouds were taken for the comparison. 

The retrieval technique, which is used in this Section is 
slightly different from that applied in Section II. In particular, 
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we derive the pair (h,/) simultaneously from the condition 
of the minimum of the cost function using reflectance spectra 
in the Og A-band (see the Appendix). We did not apply the 
technique based on the simultaneous retrieval of the pair (h, l) 
in the previous section because the retrieval method based 
on the derivation of h at the fixed value of Į gives us more 
insights in the peculiarities of the retrieval algorithm as used in 
interpretation of measurements. This was a major issue studied 
in Section II. 


Generally, we find that the agreement of results obtained 
using both instruments is quite good. In particular, both retrieval 
techniques show smaller cloud-top heights as moving from the 
equator to the north (see Fig. 7), which is a well-known effect 
related to the change in the position of the tropopause with 
latitude. There are disagreements, however (see in particular, 
Fig. 7). Let us consider obtained results in more detail, starting 
from Fig. 6. 

It follows from data given in Fig. 6 that results of the 
SACURA are inside of error bars (approximately +1 km) 
of the ATSR-2 retrievals. So, both techniques give similar 
cloud-top heights. This is an encouraging fact. 


On the same Fig. 6, we give cloud-top heights obtained 
using a very similar technique as ours but assuming that the 
cloud can be substituted by a Lambertian reflector with the 
albedo 0.8. The correspondent algorithm was developed by 
Koelemeijer et al. [9]. It is called the Fast Retrieval Scheme 
for Clouds from Oxygen A band (FRESCO). The FRESCO is 
often used to account for cloudiness in trace gas concentrations 
retrieval schemes. We see that the SACURA much improves 
the FRESCO results. This is easily understood. Indeed, we 
consider the radiative transfer inside clouds, which is not 
a primary target of the FRESCO, which can retrieve only 
so-called biased effective cloud-top heights. They are always 
smaller than actual ones (see Fig. 6). This bias is due to the 
fact that in reality photons penetrate through a cloud top and 
are absorbed by oxygen there. This leads to the increase of the 
oxygen absorption depth even for high clouds. If the process 
of the photon penetration is neglected (a Lambertian diffuser 
assumption), then this increased depth of the absorption band 
is interpreted as an existence of a cloud on a level, which lower 
that the actual altitude of a cloud. This is confirmed by Fig. 6. 
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Fig. 9. Dependence of the retrieved value of the cloud geometrical thickness 
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ATSR-2 data for the case given in Fig. 6. 


The scatter plot of cloud-top heights obtained with the 
SACURA and the FRESCO for the case given in Fig. 6 is 
shown in Fig. 8. We see that the displacement error of the 
FRESCO algorithm is 1.88 km as compared to 0.13 km for 
SACURA, if the ATSR-2 data are taken as true data. In reality, 
however, ATSR-2-derived cloud-top heights may have error up 
to 1 km and even higher. This should be remembered when 
performing comparisons as those given in Fig. 8. The root 
mean square of the cloud-top height data obtained with the 
SACURA is three times smaller than those for the FRESCO 
(0.22 against 0.66). Therefore, the example studied shows 
that the account for the cloud geometrical thickness (even as 
an effective parameter; see Rozanov and Kokhanovsky [24]) 
allows to make a reduction of the difference between retrieved 
cloud-top heights using the ATSR-2 and GOME instruments. 
Fig. 8 clearly shows limitations of the FRESCO as applied to 
cloud-top altitude retrievals. 


The SACURA also produces as an output the effective cloud 
geometrical thickness, which gives the best fit of measured and 
calculated TOA reflectance spectra in the oxygen A-band (see 
Fig. 1) for a single vertically homogeneous cloud layer. This is a 
new product, which cannot be derived from ATSR-2. However, 
as it was shown above (see Figs. 2 and 4 and also the theoretical 
study of Rozanov and Kokhanovsky [24]), the retrieved value 
of J is highly dependent on the vertical structure of the cloudi- 
ness under investigation. In particular, for multilayered cloud 
fields, l can become unrealistically large (e.g., 1 > h). This is 
the case for the orbit given in Fig. 6 as well. So, we conclude 
that most probably multilayered cloudiness has existed in the 
field of view of the GOME instrument during measurements. 
To justify this assumption we present the dependence of the re- 
trieved value of J on the degree of the coverage D of the field 
of view of GOME by cirrus clouds for the case given in Fig. 6. 
Such clouds often present near the tropopause. The values of D 
were obtained for this orbit from ATSR-2 data. The dependence 
of J on D is shown in Fig. 9 by symbols for values of D larger 
20%. It follows that the results given in Fig. 9 can be fitted by a 
linear function. This means that larger values of D lead also to 
larger values of the retrieved effective geometrical thickness /. 
This well corresponds to results of Min et al. [15]. They show 
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Fig. 10. Scatter plot of data given in Fig. 7. 


that the average photon path length increases under conditions 
of multilayered cloudiness. This results in the enhanced absorp- 
tion of photons in multilayered cloud systems as compared to 
the case of a single cloud. The main parameter, which regulates 
the absorption in our model, is J. Therefore, the increase of / 
with D is an indicator of the absorption enhancement in multi- 
layered clouds. 

We would like to underline that not for all cases we find such 
a good agreement of the SACURA and the ATSR-2 retrieval 
results as in Fig. 6. The example is shown in Fig. 7, where com- 
paratively large deviations are found (e.g., see pixels close to 
equator and 60°S latitude). Generally, deviations decrease for 
cases with | < h. For instance, the deviations Ç for pixels lo- 
cated in the range of latitudes 3.2°N to 45.3°S (where l < h) 
are equal to 0.5 km in average as compared to ¢ = 1.0 km for 
pixels located in the range of latitudes 45.3°S to 64.2°S (where 
l >h). 

The scatter plot of data given in Fig. 7 is shown in Fig. 10. It 
follows from Fig. 10 that the largest discrepancy takes place for 
low and high clouds. 

However, we suspect that although all pixels in Fig. 7 are 
completely covered by cloudiness, main assumptions behind ei- 
ther the SACURA or the ATSR-2 retrieval technique are vio- 
lated. This may explain the discrepancy found. 

In particular, let us consider pixels in the vicinity of the 
equator (see Fig. 7), where the SACURA gives too low values 
of h. This can be explained by the fact that cloud altitudes 
are highly variable near the equator. So, the assumption of the 
horizontally uniform layer, used by the SACURA is violated 
in this area for large GOME pixels (even for the case of 100% 
cloudy pixels). Then, measurements of h with the GOME 
instrument become inadequate. Our guess is supported by the 
ATSR-2 data given in Fig. 7 (the points close to right vertical 
axis), which show a high variability in a cloud-top height on a 
small spatial scale. We also retrieved the cloud optical thickness 
for this area using algorithm described by Kokhanovsky et 
al. [12]. The value of the cloud optical thickness varies in the 
range 17—107 on the distance of approximately 80 km. This 
confirms our guess that the cloud field considered is highly 
inhomogeneous in the horizontal direction. Therefore, we face 
here problems not related to the SACURA itself but to the 
GOME spatial resolution. 
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Fig. 11. Dependence of retrieved cloud-top height on latitude obtained using 
different retrieval techniques and optical instruments onboard the ESR-2. 
SACURA/GOME results are given by the solid line with crosses. Parallax and 
infrared ATSR-2 techniques are shown by symbols. The dashed line shows the 
dependence of the geometrical thickness excess (1 — h) on the latitude. Results 
for the infrared ATSR-2 technique and the SACURA/GOME retrieval coincide 
with those given in Fig. 7 around 60°S. 


Let us consider now the disagreement of retrievals close to the 
60 °S latitude (see Fig. 7). The discrepancy between values of h 
obtained from the GOME and the ATSR-2 is around 2 km in this 
region. Whether this is due to the failure of the SACURA, the 
ATSR-? retrieval technique or both remains unclear. However, 
we would like to underline that generally the accuracy of the 
infrared ATSR-2 cloud-top height retrieval technique decreases 
for low-level clouds (especially if cirrus clouds are present in 
the scene). Then, the error of the infrared technique can reach 2 
km [19], which can explain the discrepancy found. 

The comparison between GOME SACURA and ATSR-2 in- 
frared and parallax [21] techniques is given in Fig. 11 for the re- 
gion ~ 60 S. Note that data for SACURA/GOME and infrared 
ATSR-? technique coincide with those given in Fig. 7. The par- 
allax technique uses the dual-viewing capability for the same 
scene. It gives values of h that differ from those obtained both 
from cloud brightness temperature measurements performed by 
ATSR-2 and from SACURA results for this particular case. 

Let us underline that the retrieved effective geometrical thick- 
ness / is larger than h for the cases studied. To illustrate this fact, 
we show the geometrical thickness excess € = l — h in Fig. 11 
for latitudes 57° to 65°S. Note that € is positive and reaches 3.5 
km at 64°S. This points again to the fact that the multilayered 
cloudiness probably existed in the field of view of the detector, 
and the assumption of a single homogeneous cloud layer used 
in the retrieval procedure is most probably wrong. 


IV. CONCLUSION 


We presented here the first application of the SACURA to the 
determination of the cloud-top height h and geometrical thick- 
ness / using the TOA reflectance measurements in the oxygen 
A-band. The algorithm is relied upon the semianalytical solu- 
tion of the radiative transfer equation, which is valid for single 
homogeneous clouds having optical thickness larger than 5. 

Using the experimental and simulated reflectance spectra in 
the Op A-band, we have shown quantitatively the dependence of 
the retrieved value of A on the geometrical thickness of clouds 
l. We have demonstrated that parameters h and Į can be re- 
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trieved from the measured reflectance spectra in the oxygen ab- 
sorption band simultaneously (in the assumption of a homoge- 
neous layer). This is similar to findings reported by Daniel et al. 
[3], who used, however, not only oxygen A-band but also other 
oxygen absorption bands (along with absorption bands of the 
O2—-Op collision complex). 

We have found unrealistically large values of the retrieved 
geometrical thickness (l > h) for most cases studied. Our the- 
oretical study shows that such large retrieved values of J under- 
line the fact that a multilayered cloudiness existed in the field of 
view of the instrument during measurements. 

Unfortunately, the direct measurement of the vertical struc- 
ture of the cloudiness studied is not available to us. Neverthe- 
less, the found relationship between retrieved values of | and 
the degree of the coverage of the field of view of the GOME by 
cirrus clouds points out that our assumption is quite reasonable. 

Although results obtained as far as the value of / is concerned 
are of a limited value, we believe that we can identify the cases 
with multilayered cloudiness. 

The new information which can be derived from the re- 
flectance spectra in the oxygen A-band is not very significant 
in the case of such large ground pixels as that by GOME 
(40 x 80 km? in the best case). Then, the multilayered structure 
of cloudiness is quite probable to occur. But our technique could 
be of importance for future spaceborne optical instruments 
including the Orbing Carbon Observatory. 

The large retrieved values of | enable us to reach coinci- 
dence between the measured and simulated A-band spectra of 
the real-world multilayered cloud systems and a single homo- 
geneous layer model. Therefore, the retrieved cloud-top height 
is more realistic comparing with other techniques, which com- 
pletely neglect multiple light scattering in the cloud itself or ac- 
count for multiple light scattering using fixed value of the cloud 
geometrical thickness. 

The comparison of the obtained cloud-top height with other 
results as retrieved from infrared ATSR-2 and METEOSAT 
measurements shows that our algorithm is capable to give 
the cloud-top height similar to these popular techniques. The 
discrepancy is smaller than 1 km for the most cases. 


APPENDIX 
INVERSE PROBLEM SOLUTION 


The theoretical basis of the SACURA is given by Rozanov 
and Kokhanovsky [24]. However, for the sake of clarity and 
completeness, we summarize the main steps of the inverse 
problem solution here. 


A. Linearization of the Reflectance 


Fist of all, the TOA reflectance R is presented in the form of 
a Taylor expansion around the assumed value of the cloud-top 
height equal to ho 


R(h) = R(ho) + S ai(h — ho) (1) 


where a; = R (ho) /i!. Here R® (ho) is the i-derivative of R 
at the point ho. The next step is the linearization, which is a stan- 
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dard technique in the inversion procedures [22], [23]. We found 
that the function R(h) is close to a linear one in a broad interval 
of the argument change [11]. Therefore, we neglect nonlinear 
terms in (1). Then, it follows that 


R(h) = R(ho) + R'(ho)(h — ho) (2) 


where R’ = (dR/dh). We assume that Rè is measured at several 
wavelengths in the oxygen A-band. Then, instead of the scalar 
quantity R, we can introduce the vector Ř nes with components 
(R(à1), R(à2), - . ., R(An)). The same applies to other scalars 
in (1). 

Therefore, (2) can be written in the following vector form: 


y= ax (3) 


where 7 = Rimes — R(ho). q= R'(ho), and z = h — hg. Note 
that both measurement and model errors are contained in (3). 


B. One-Parameter Retrieval Algorithm 


The solution ĉ of the inverse problem is obtained by the min- 
imizing the following cost function [22]: 


® = ||7 - x|]? (4) 


where || || means the norm in the Euclidean space of the corre- 
spondent dimension. 

The value of ĉ, where the function ® has a minimum can be 
presented as 
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where (,) denotes a scalar product in the Euclidean space, and n 
is the number of wavelengths, at which the reflection function 
is measured. 

Therefore, knowing the measured spectral reflection func- 
tion Rimes and also the calculated reflection function R(ho) and 
its derivative F’(ho) at several wavelengths, the value of the 
cloud-top height can be found from (5) and equality A = ĉ+ ho. 
The value of hg can be taken equal to 1.0 km, which is a typical 
value for low-level clouds. 

The reflection function of the cloudy atmosphere depends 
not only on the value of h but also on / in the oxygen absorp- 
tion band. In particular, / is a major parameter that determines 
the level of light absorption in cloud in the oxygen absorp- 
tion band. Therefore, the value of ê depends on the assumed 
cloud geometrical thickness J;. Accordingly, the cost function 
P(A(L)) (AC) = ê(l:) + ho) depends on the choice of l; as 
well. 

The cost function at a given argument l; characterizes the de- 
gree of the coincidence of the model and the experiment for the 
derived value of ĉ(l;). Therefore, it is natural to take as retrieved 
values of the cloud geometrical thickness and cloud-top height 
those values for which ®(h(J;)) reaches a minimum. 

The algorithm of the determination of values h and I is re- 
duced, therefore, to the choice of l; (i = 1, . . . , N), the estima- 
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tion of ĉ(l;) for each l; in correspondence with (5) and the de- 
termination the value of ¿, which insures that ®(A(l;)) — min. 

The technique described above is used for the solution of the 
problem under consideration in Section II. 


C. General Retrieval Technique 


In Section HI, both values of h and } are found simultane- 
ously. This requires the minimization of the following cost func- 
tion [see (4)]: 


ò = ||g— ÂX. (6) 


The elements of the matrix A are correspondent weighting func- 
tions [22]. The solution of the inverse problem is given by the 
vector-parameter X. This vector has five components in our 
case, which give corrections to the initially assumed cloud-top 
height and cloud geometrical thickness, the correction to the ini- 
tially assumed half-width of the spectrometer spectral response 
function, the shift parameter, and the squeeze parameter. 

Clearly, the first two parameters give us final values of the pair 
(h, 1) to be retrieved. The third parameter allows to adjust the 
assumption on the instrument response function. The last two 
parameters allow the reduce errors related to the displacement 
of the experimentally measured spectrum due to the errors of 
the spectral calibration. 

We have developed two versions of the retrieval algorithm. 
One is based on the exact radiative transfer calculations of the 
reflection function R, and another one is based on the approxi- 
mate representation of R by the following equation: 


R= Ro + DRT (7) 


where Ro gives the reflection function of the part of atmos- 
phere above the cloud, Re is the cloud reflection, and functions 
T; (i = 1,2) give transmission coefficients from the sun to a 
cloud and from the cloud to a satellite, respectively. Approxi- 
mate equations for all functions in (7) are given elsewhere [11]. 
Note that we have used the retrieval technique based on the ap- 
proximate representation of R [see (7)] in this paper. It allows 
to speed up the retrieval process considerably. 
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